Generalized theory for node disruption in finite size complex networks 
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After a failure or attack the structure of a complex network changes due to node removal. Here, 
we show that the degree distribution of the distorted network, under any node disturbances, can be 
easily computed through a simple formula. Based on this expression, we derive a general condition 
for the stability of non-correlated finite complex networks under any arbitrary attack. We apply 
this formalism to derive an expression for the percolation threshold fc under a general attack of the 
form fk ^ k'^ , where fk stands for the probability of a node of degree k of being removed during the 
attack. We show that fc of a finite network of size A'' exhibits an additive correction which scales 
as A''"^ with respect to the classical result for infinite networks. 

PACS numbers: 89.75.Hc,89.75.Fb,02.50.Cw 



The stability of graphs against various disrupting 
events is a central issue in the study of complex net- 
works [H, S, S, i, i, S, 0, S, i, E^. If information is 
transported across a network, as is the case of epidemics 
across social networks or information broadcast through 
Internet, the "damage" of some nodes can dramatically 
affect the dynamics of the system. In the context of dis- 
ease spreading, this could lead to an epidemic extinction, 
while in communications to a halt of information broad- 
cast 11, 12, 13, 14, 15[. Both, the topological structure 
of the network and the nature of the attack determine the 
resulting effect [l^. For example, it has been shown that 
scale- free (SF) networks display a high degree of tolerance 
against random failures [9|, while, on the other hand, 
they are quite sensitive against intentional attacks [lo| . 
Clearly, there are various strategies to perform an inten- 
tional attack, and each one of them requires a different 
level of knowledge on the network top ology [IgI . 17. l^. 
A rather general attack, proposed in [1^, |20| and that 
we will use in this paper, takes the form ^ k'^ , where 
fk denotes the probability of a node of degree k of be- 
ing removed during the attack, while 7 is associated to 
the degree of knowledge of the attacker. The analysis of 
this attack has revealed that in SF networks an increase 
of 7 leads to a decrease of the critical fraction of nodes 
that must be removed to disintegrate the network, i.e., a 
decrease in the percolation threshold fc H, 



Though many results have been derived for infinite SF 
networks, very little is known about the stability of finite 
networks. Typical examples of small size finite networks 
are ad-hoc networks of commercial mobile devices, fre- 
quently used for communication [21|] , temporary peer-to- 
peer networks formed by Bit Torrent clients for efficient 



download of file [24I and networks of autonomous mo- 
bile robots [ill . The operation of these systems relies on 
the robustness of the highly dynamical underlying net- 
work. Thus, a good understanding on the stability of 
these small size networks is imperative for these applica- 
tions. Moreover, we can say that in general a comprehen- 
sive theory for the stability of arbitrary finite networks 
under any node disturbance is still lacking. 

In this paper we attempt to shed some light on this 
matter by proposing an alternative derivation for the 
percolation threshold. Let us clarify that in this paper 
we are interested in understanding only the emergence 
of a giant component at percolation threshold. Hence, 
the term percolation and emergence of giant component 
can be considered synonymous. In our approach, in- 
stead of applying a generating function formalism to find 
an analytic expression for the percolation threshold as 
in [ll,[2,[l3|, we used the fact that during an attack the de- 
gree distribution of the network changes (see Fig.[T]). We 
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FIG. 1: The scheme illustrates an attack as consisting of two 
steps: (a) selection of nodes to be removed (see set R), and 
(b) cutting of the edges that run from the surviving nodes 
(represented by set S) to the set of removed nodes. As the 
scheme shows, the attack affects also the degree of the sur- 
viving nodes. 
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show that the degree distribution of the distorted (un- 
correlated) network, under any node disturbances, can 
be easily computed through a simple formula. Based 
on this expression, we derive a general condition for the 
stability of non-correlated complex networks under any 
arbitrary attack. This condition applied to the study 
of network stability under the general attack proposed 
in 1^, 23| leads us to a general expression for the perco- 
lation threshold fc- We show that fc of a finite network 
of size N exhibits an additive correction which scales as 
with respect to the classical result for infinite net- 
works [l|, 0, [l^ . Simulation results confirm all these find- 
ings. 

Network topology after disturbance. A failure or attack 
can be thought in the following way. Let pk be the degree 
distribution of the network before the attack. The first 
step in the attack is to select the nodes that are going 
to be removed. Let us assume that this is performed by 
means of fk, where fk represents the probability for a 
node of degree k being removed from the network. Note 
that the only restriction on is < < 1. After 
the node selection, we divide the network into two sub- 
sets, one subset contains the surviving nodes (S) while 
the other subset comprises the nodes that are going to 
be removed (R). Fig. [1] illustrates this procedure. At 
the moment the nodes in R are actually removed, the 
degree distribution of the S-nodes is changed due to the 
removal of the E edges that run between these two sub- 
sets. The probability (j) of finding an edge in subset S 
that is connected to a node in subset R is expressed as: 



X]i=0 ^ Pi fi 



{j:Z,kpk)~l/N 



(1) 



The reasoning behind this expression is as follows. The 
total number of half-edges in the surviving subset, in- 
cluding the E links that are going to be removed, 
is Yl,%.oj{^Pj){^~fj)- The probability for a ran- 
domly chosen half-edge of being removed is simply 

Ej^o * i^Pi) /^/(Efelo k i^Pk) - !)■ E is the number 
of half-edges in S times this probability, and (j) is obtained 
by dividing E by the number of half-edges in the subset 
S. Notice that the removal of nodes can only lead to a 
decrease of the degree of a node. Finally, to calculate the 
degree distribution pj. after the attack, we still need to 
estimate the probability of finding a nodes with degree 
q in the surviving subset S (before cutting the E edges). 
This fraction takes the simple form: 



^~YT=0Pif^ 



(2) 



Now we are in condition to compute pj,. Using Eqs. ([T]) 
and ([2]), we obtain the following expression for p'f.: 



p'k = E 

q—k 



iq-k 



(3) 
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FIG. 2: Degree distribution before (circles), and after (stars) 
the attack. In (a) the network topology corresponds to Erdos- 
Renyi graphs, with (fc) = 5, while in (b) to scale-free networks, 
Pk ~ fc~° with a = 2.5 [23|. Symbols correspond to simula- 
tions, solid lines to the theoretical p^, given by Eq.(|3}, and 
dashed curves to the (theoretical) initial pk- In (a) and (b) 
the main figure corresponds to networks with TV = 10^ nodes, 
while the inset shows the result for small networks with just 
iV = 50 nodes. 



Eq. ([3]) can be iteratively evaluated by replacing pk with 
p'f. into Eqs. ([T]) to ([3]). It is instructive to notice that for 
failure, i.e. fk = /, and assuming N :s> 1, Eq. ([T]) reduces 
to (j) = f while Eq. ^ becomes p^ — Pq- In consequence, 
from Eq. ([3]) we retrieve the degree distribution pj, after 

failure which reads p'k = Yl^=k (y\ 

A similar expression has also been used to described pj, 
after an ad-hoc attack in scale-free networks with N » 

I M. 
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Fig. [5] shows a comparison between stochastic simula- 
tions (symbols) and Eq. dS]) (solid line) for two different 
network topologies, namely Erdos-Renyi (ER) graph (a) 
and scale-free (SF) networks (b) [25|]. Removal of nodes 
(and edges) was performed through an attack of the form 
fk ~ A:''', with 7 = 1. In the figure two different system 
sizes are shown: TV = 10^ and TV = 50 (figure insets). 
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Critical condition. The following expression tells us 
whether an infinite network percolates after an attack Q : 



(k)' 



> 2. 



(4) 



where {k)' and (fc^)' refer to the first and second mo- 
ments of the degree distribution after the attack. We 
borrow the critical condition for infinite networks given 
by Eq. ([4]) to define a "percolation" criterion for finite 
networks. Thus, by definition we assume that the con- 
dition kI — 1 determines the point at which the network 
breaks down [Til, l20l|. To compute (fc)' and (fc^)', we 
utilize the generating function: 
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(5) 



After exchanging the order of the sum, the Binomial the- 
orem can be applied, and we obtain: 



(6) 



FIG. 3: Percolation threshold fc under changes of the at- 
tack exponent 7 for three different SF networks, pk ~ 
with a = 2,2.5 and 3, and A'^ = 10^. Symbols correspond 
to stochastic simulations, while solid curves correspond to 
Eqs. ^ and (|lUp . The horizontal dashed lines indicate the 
asymptotic value of fc given by Eq. (|lip . 
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From Eq. ([6]), the first two moments can be easily com- 
puted as (fc)' = dGo{l)/dx and {P)' = d'^Go{l)/dx^ + 
dGo{l)/dx. After some algebra we obtain that the criti- 
cal condition given by Eq. (|4]) takes the form: 

\Tpk{l - fk)k] n^Pkil - fk)k' +Y,Pk{fk - 2)fc j 

+ (^Pk{l-fk){2-k)k^ =0(7) 

Eq. ([7]) determines the stability condition (according to 
the given definition) for any uncorrelated network of fi- 
nite size under any arbitrary attack. In the limit of 
N ^ 00, Eq. ([7]) reduces to: 



fc of nodes that is required to remove through such at- 
tack in order to destroy the network, i.e., the percolation 
threshold. Thus, the problem reduces to compute, for 
a given 7, the critical attack intensity C* . Replacing in 
Eq. ([7]) the above definition of fk and after some algebra 
we obtain: 

^* = 2^W^) ^ ['^^^'^ ^'/^^ ^'^ + 

1/2 



+ H/N) {2Q +{q+ (fc^+i))') - 4gp(fc^+i) 



where p = (fc^) - 2(fc), q ^ (fc^+i) - {k'<+^), Q = 
(k'^+'^)p+ (k)q, and (fc") is defined as (fc") = J2k k"" Pk- 
Since the fraction of removed nodes / after an attack is 
/ — '^j.Pkfk, the expression for percolation threshold fc 
is simply: 



Y,Pkk{k{l- h) + fk-2) = Q. 



(8) 



/c = C*(fc^). 



(10) 
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Interestingly, Eq. ^ can be also derived through a more 
classical generating function formalism [il ]. 

Generalized attack. In the following, we model the vari- 
ous dynamics (attacks) through a generalized equation of 
the form: fk = C k^ , where 7 is a real number signifying 
the amount of network structure information available 
to the attacker to breakdown the network [2^, and G is 
a constant that we refer to as attack intensity. Clearly, 
7 > represents a situation in which high degree nodes 
are removed with higher probability, while 7 < models 
the opposite. The last case is suitable to situations in 
which low degree nodes are more prone to fail. We are 
interested in knowing, for a given 7, the critical fraction 



Fig. [3] illustrates the behavior of fc on three SF networks 
[a — 2, 2.5 and 3) upon changes in the attack exponent 
7. The symbols correspond to stochastic simulations per- 
formed on the networks of size N — 10^, while the black 
curves refer to Eqs.® and (fTO|) . In the numerical ex- 
periments we have computed fc following (20j : when the 
fraction of removed nodes is /c, the probability F of find- 
ing the network with k' > 2 is 1/2. 

It is interesting to observe that for any SF network, 
the minimum fraction <I>c of nodes that is required to be 
removed to break down network is obtained by taking 
the limit 7 ^ 00 of Eq. 



(11) 



$c = lim fc{l,a) = h{a)- — jr^ 
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FIG. 4: Percolation threshold fc as function of A^. Symbols 
correspond to stochastic simulation, solid curve to Eq. (I12|) . 
and dashed curve to the classical result f^. The error bars 
represent the confidence interval that contains f^{N). Simu- 
lations were performed with ER networks with average degree 
(k) — 3 and maximum degree fci/ ~ 5. The inset shows the 
scaling of — fc{N) with respect to A''. The best fitting of 
the data corresponds to a slope —0.92 ± 0.05 (dashed line). 



where h{a) is h{a) — {k"^) ~2{k) and /cm is the maximum 
degree of the original network. Notice that Eq. (fTT|) rep- 
resents an attack performed having full knowledge on the 
network topology. The asymptotic values corresponding 
to Eq. ([TT|) are shown in Fig. [3] as horizontal dashed lines. 
Fig. [3] indicates that typically an increase of the infor- 
mation about the network topology, resp. 7, helps the 
attacker to break down the network with the removal of a 
smaller number of nodes. However, information becomes 
redundant as the asymptotic value <i>c is approached. 

Finite network size. To illustrate the effect of network 
size N upon the percolation threshold fc{N), we cus- 
tomize Eq.® for random attack or failure. When 7 = 
and fk is independent of fc, we obtain: 



UN) = + 



1 

N 



-m/{k)) \ 
{k^)-{k) J 



(12) 



where is the well-known percolation threshold for in- 
finite networks under failure P, i] which reads: = 
l-[l/((fc2>/(fc}-l)]. 

Fig |4] shows a comparison between Eq. (|12p (solid 
line) , (dashed line) , and stochastic simulations (sym- 
bols) for ER networks of different sizes. Notice that 
Eq. IT^ predicts the correct scaling of fc with N, i.e., 
fc{N) — r-^ N^^. The observed deviation between 
Eq. p2| and simulations can be arguably attributed to 
correlations effects, which have been ignored in the cur- 
rent approach. 

Concluding remarks. We have proposed a general pro- 
cedure to calculate the distorted degree distribution of 



uncorrelated finite size networks under arbitrary fail- 
ure/attack. Using the expression for the distorted degree 
distribution we have derived the critical condition for the 
stability of finite size networks. The formalism has been 
further applied to derive an expression for the percola- 
tion threshold under a general attack. Finally, it was 
shown that the obtained percolation threshold predicts 
an additive correction which scales as with respect 
to the classical result for infinite networks, as observed 
in simulations. 

The results derived throughout this manuscript are 
valid only for uncorrelated networks. The effect of cor- 
relations on the percolation threshold for finite (and in- 
finite) networks remains as one of the major challenges. 
Further extensions of this theory will be focused in that 
direction. 
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